A 2-D staggered grid pseudospectral and finite difference hybrid method is applied to perform numerical simulations of lunar seismic wave propagation in a laterally heterogeneous crustal Moon model. The relief of lunar crust and mantle interface is defined with different height, width and position based on the present lunar structural models. The effects of these factors on lunar records are discussed, respectively. How the lateral variations of Moon crustal thickness affects the strong wave coda of lunar signals is revealed. The results show that the height of the relief dominates the formation of strong wave coda, greatly affects the duration time and amplitude of the coda. The width of the relief has impact on the amplitude. The effect of position of the relief is complex, not only affecting the travel time of P and S waves, but also the duration time and amplitude of the coda.